Exactly solvable two-level quantum systems and Landau-Zener interferometry 
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I present a simple algorithm based on a type of partial reverse-engineering that generates an unlim- 
ited number of exact analytical solutions to the Schrodinger equation for a general time-dependent 
two-level Hamiltonian. I demonstrate this method by deriving new exact solutions corresponding to 
fast control pulses that contain arbitrarily many tunable parameters. It is shown that the formal- 
ism is naturally suited to generating analytical control protocols that perform precise non-adiabatic 
rapid passage and Landau-Zener interferometry near the quantum speed limit. A general, exact 
formula for Landau-Zener interference patterns is derived. 



Although they have pervaded quantum physics since 
its inception, very few time-dependent two-level quantum 
systems are known to be analytically solvable. Among 
the most famous examples of exactly soluble two-level 
evolution is the Landau-Zener-Stuckelberg-Majorana 1 " 5 
(LZSM) problem, which remains a very active area of re- 
search due to numerous applications pertaining to quan- 
tum phase transitions^ quantum control 7- — and quan- 



tum state preparation. 



The hyperbolic secant pulse 



of Rosen and Zener comprises another well known ex- 
act solution discovered in the same period^ This so- 
lution has played an important role in the contexts of 
self-induced transparency's and qubit control^ 7 - - — and a 
number of authors have since found that this example be- 
longs to a larger family of analytical controls.— ~— In the 
last few years, several of these examples have proven very 
beneficial to the fields of quantum control and quantum 
computation ) 17 ! 18 ' 31 " — where analytical solutions are of- 
ten central in the design of control fields that are fast, 
precise, and robust against noise. However, the rarity 
of exactly solvable two-state systems has severely lim- 
ited one's options in developing an analytical approach 
to qubit gate design. 

In a recent work,— a systematic method for deriving 
arbitrarily many families of exactly solvable two-state 
systems was presented, vastly extending the number of 
known analytical solutions. This method allows one to 
input many of the basic features of the desired control 
field and then compute exactly the corresponding evolu- 
tion of the system with the provided formulas. However, 
a limitation of this work is that it applies only to sys- 
tems where the driving is along a single axis of the Bloch 
sphere, such as in the case of electrically-driven singlet- 
triplet qubits,— ~— making it inapplicable to the majority 
of driven two- level systems. 

In this paper, I address this limitation by presenting a 
method to generate arbitrarily many families of solutions 
in the most general case where the two-level Hamiltonian 
has time-dependence along any set of axes of the Bloch 
sphere. Of course, one can easily generate exactly solv- 
able Hamiltonians by first choosing the evolution opera- 
tor and then differentiating to obtain the corresponding 
Hamiltonian, but it is challenging to arrive at a phys- 
ically meaningful Hamiltonian in this way. In contrast, 



the method presented here allows one to specify the basic 
form and many of the features of the Hamiltonian whose 
evolution one wishes to solve before proceeding to com- 
pute the exact solution for this evolution. This method 
for generating exact solutions has important applications 
in a vast range of quantum physics problems, includ- 
ing the development of quantum controls for essentially 
any quantum computing platform and control protocols 
for performing LZSM interferometry and non-adiabatic 
rapid passage (NARP). I illustrate this by deriving new, 
exactly solvable LZSM driving fields and control pulses 
that execute a desired evolution at speeds approaching 
the quantum speed limil42r— (QSL). Attaining fast evo- 
lution times is especially crucial in quantum computing 
where quantum gates need to be performed on timescales 
much shorter than the decoherence time. In the case of 
periodic driving through a level anti-crossing, I show that 
the formalism allows one to easily derive analytical ex- 
pressions for LZSM interference patterns and conditions 
for coherent destruction of tunneling ) 45 ' 46 

The Hamiltonian we will consider has the general form 



H = b(t) ■ a, 
b(t) = [b x (t),b y (t),b z (t)} 



(1) 



where b(t) = [b x (t),b y (t),b z (t)] is a three-vector whose 
components are arbitrary real functions, and a = 
[fx j <?Vj o'z] is a vector of Pauli matrices. This Hamil- 
tonian describes any time-dependent two-level system. 
The components of the effective magnetic field b(t) can 
each admit various physical interpretations, most com- 
monly as cither driving fields or time-dependent energy 
splittings. The evolution operator for this Hamiltonian 
obeys a Schrodinger equation whose form can be made 
compact by first parametrizing the evolution operator as 
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\u 2 i\ 2 = l, (2) 



and then transforming to a rotating frame: v±i = 
e if*dt'b z (t') Uii ^ V2i = e -iftdt'b z (t') U2 ^ Writing i x = 

/?(<) cos <p(t) and b y = j3(t) sin tp(t), and defining a(t) = 
2 J* dt'b z (t') - tp(t), we have 



i'21 



(3) 



Here, it is manifest that the evolution operator in the ro- 
tating frame depends on only two real functions, a and 



2 



/?. We must further specify b z to return to the lab frame, 
but this choice can be made after the evolution opera- 
tor is computed in the rotating frame. In what follows, 
I will show that one can systematically find analytical 
solutions with either a or j3 chosen as desired; although 
one cannot completely specify both a and /3 at will (were 
this the case, one could solve all possible two-state prob- 
lems analytically), we will see that one still has a large 
amount of control over the features of the second, un- 
specified function. 

For concrctcness, we will suppose that one wishes to 
fix P(t) at the outset (the formalism can easily be mod- 
ified to fix a(t) instead). While we cannot then find 
analytical solutions for arbitrary a, there exists a dif- 
ferent parametrization of the Hamiltonian in which a is 
replaced with a new function, Ki(t), such that one can 
systematically generate an analytical expression for the 
evolution operator for arbitrary choices of /3 and kj. To 
see this, first express the rotating- frame evolution oper- 
ator in terms of some complex K,(t): 



1 ° ) 



V21 = -irj e 



-iJ*dt'P(t')e-< t "> 



(4) 

with r\ = ±1. This choice of parametrization is motivated 
by observing that we can combine the two equations in 
([3]) to obtain /3 2 = — (vn/vn)(v2i/v2i), which implies 
i>n/vn = —ifie K , 1)21/^21 = —i[3e~ K for some complex k. 
Consistency between Eqs. ^ and ([3]) requires 



a(t) = -in(t) +r]- -2 



dt'P(t')smhn(t'), (5) 



which should be interpreted as follows: For any choice of 
a complex n(t) and real /3(t) such that the a(t) computed 
from Eq. ([5]) is real, the evolution operator obtained from 
Eq. (j4} is the exact solution for this a and (3. Writing 
k = Kfj + iki and imposing Im(a) = determines kr in 
terms of kj: kr = — 2 tanh -1 tan (x + 7r/4), with 



dt'p(t')sm(Ki(t')). 



(6) 



From this we find an expression for a that is real for any 
choice of ki: 



a = m + ri--2 dt'P(t') cos K^t') cot [2x(t')}. (7) 
2 Jo 

While this parametrization has the nice feature that 
kj can be chosen freely, one drawback is that one must 
then perform the integration in Eq. ([6]) , making it harder 
to relate the features of kj to the driving field a. We can 
avoid this by specifying x directly, but at the expense of 
now having to choose functions \(t) that obey |x| < \P\, 
which arises directly from Eq. ((6|). In terms of x-, the 
lab-frame evolution operator and driving field a are 
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(8) 
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2/Vl-x7/3 2 cot(2 x ). (9) 



The initial conditions Uii(O) = 1, W2i(0) = imply 
X(0) = 0, and x(0) = —77/3(0) ensures d(Q) is finite. 
Eqs. dSJ) , (J9j) embody one of the main results of this pa- 
per, as they constitute a general analytic solution of the 
evolution generated by the Hamiltonian of Eq. ((T|) . The 
task of finding analytical solutions has been reduced to 
first picking j3(t) at will. One then selects x(t) to pro- 
duce a desired a(t) via Eq. ©, fixing the rotating- frame 
Hamiltonian. Once these choices are made, an analytical 
expression for the evolution operator follows immediately 
from Eq. {8]). The fact that x directly controls \uu\, 
\u2i\ makes this formalism a powerful tool in designing 
analytical driving fields that execute precise quantum 
operations or LZSM interferometry, as will be demon- 
strated shortly with explicit examples. A simple case is 
X(i) = -r)f*dtfp(tf), where Eq. © gives that a = 0, 
corresponding to an x- rotation for any /3. Another sim- 
ple example arises when (3 = x — 0, for which Eq. ([8]) 
yields a z-rotation in the lab frame for any choice of b z . 

The physical origin of the constraint |x| < \/3\ lies in 
the notion of the quantum speed limit ^t£l which refers 
to the minimum time it takes a quantum state to evolve 
to a different state in the Hilbert space due to energy- 
time uncertainty. Indeed, |x| < \j3\ implies that the 
fastest possible evolution from x(0) = to x(T) > 
is obtained by choosing x(^) = f dt'\/3(t')\, with the 
shortest time given by substituting t = T in this ex- 
pression and solving for T. For constant f3 — (3q > 0, 
we immediately obtain T m i n = x(^ 1 )//3o J which is the 
QSL time Tqsl derived for an arbitrary Hilbert space 
in Ref. 43. We will refer to |x| < as the QSL con- 
straint. The present work leads to a general definition 
of T QSLl X {Tqsl) = Jo QSL dt'\/3(t% for arbitrary time- 
dependent two-level systems. Notice that the QSL evo- 
lution x = ± Jq dt' j3(t') coincides with a = 0, suggesting 
that the fastest quantum operations are those which tend 
to minimize d, a tendency that is borne out in the exam- 
ples given below. 

The fact that the QSL appears as a simple condition 
on x makes the formalism of Eqs. JS]),© very effective for 
designing quantum controls that operate near the QSL. 
To see how this works for a general Pit), note that a 
simple way to construct a function x(^) which obeys the 
QSL constraint is to first find a function which satis- 
fies the constraint in the case where (3(t) = fio is a con- 
stant. Denoting this latter function by Xo(^) and defin- 
ing B(t) = J* dt' (3(f), if we choose x (t) = Xo{B(t)/p ), 
then |x| = \Pxo(B/Po)/Po\ < \P\ automatically follows. 
[Note that all the single-axis driving examples of Ref. l35l . 
where the notation there is related to the present by 
q = cos(2x), can be extended to multi-axis solutions 
using this trick.] Furthermore, if the control field cor- 
responding to xo operates near the QSL, this will also 
tend to be the case for the one generated by x- Focusing 
then on the case pit) = Po,we can construct controls that 
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FIG. 1: (Color online) Control field b z generated by the x 
from (a) Eq. ([I0j) with b x = /3o, ip = 0, k = 6, a 2 = 0, 
a 6 = 4/(ttA)), (b) Eq. (H) with b x = #,(1 + sin 2 (2/3 t)/2), 
V? = 0, fc = 6, ci2 = 04 = 0, ae = 4/7T. A Hadamard gate is 
achieved for total evolution duration T = 3.61/ /?o- 



operate near the QSL by choosing a x(i) which contains 
parameters that can be tuned to values where x = ^fiot 
and the QSL constraint is saturated. 

To illustrate this method, consider the example 

x = + M) 2 + (a^t) 4 + ■■■ + M) fc r 1/fe , (io) 

where /3(t) = /3q and the a% are arbitrary constants, A; 
is an even integer, and rj = 1. The QSL constraint is 
satisfied regardless of how large k is, so that this x yields 
an exact solution with arbitrarily many parameters dj. 
We can make the corresponding control field a pulse by 
setting a-k — 4:f3 /ir, so that x ~> ~tt/4 and a — > as 
t — > oo. 52 Examples of these pulses are shown in Fig. [TJi. 
The duration of the pulse approaches Tqsl in the limit 
Oi<fe - > 0, fc — > oo, as can be seen by observing that 
X —> —0ot in this limit. The substantial amount of tun- 
ability in this solution already makes it very attractive 
for applications in quantum computation such as dynam- 
ically corrected gates) 32 ' 34 ' 39 where one tunes the shape of 
the pulse to perform a specific quantum operation while 
simultaneously suppressing leading-order errors. 

Using the prescription outlined above, we can extend 
this solution to the case of non-constant f3: 

X = -B{t)[l+{a 2 B{t)) 2 + {a i B{t)f+- ■ ■+(a k B(t)) k ]-^ k . 

(11) 

This class of pulses can be used to implement quantum 
operations by tuning b z (t) while keeping the functions 
b x (t) and b y (t) fixed. We will demonstrate this by design- 
ing a fast b z pulse that, together with b x , implements a 
Hadamard gate, a quantum operation that is ubiquitous 
in the field of quantum information processing and which 
is equivalent to a 7r-rotation about x + z. First choose 
a/j = 4/tt, which ensures that |un|, |ti2i| — > l/v2- Sup- 
posing ip — 0, if we let the system evolve for a time T such 
that Jq dt'y/ j3 2 — x 2 csc(2x) = —5ir/4, then the phases 
of un and U21 will also attain their Hadamard values. An 
example of one such b z pulse is shown in Fig. []J> for the 
case where b x is a particular oscillatory function. From 
the figure, it is evident that the b z pulse quickly sets the 
magnitudes of u\\ and u^x, while the remainder of the 
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FIG. 2: (Color online) Diabatic and adiabatic energy levels 
as a function of b z . 

evolution with b z ~ allows their phases to accumulate. 
As before, the duration of the pulse approaches Tqsl as 
Oj<fe — )• 0, fc — > oo. This example illustrates how one can 
use this formalism to design analytical single or multi- 
axis quantum controls near the QSL in the presence of 
additional continuous-driving fields. 

In addition to designing control pulses that imple- 
ment fast quantum gates, the present formalism is also 
natural for designing driving fields that perform con- 
trolled LZSM interfcrometry and non-adiabatic rapid 
passage^ 3 phenomena which have many applications in 
quantum control^— state preparatio n 12 ' 13 and qubit 
readout^- 3 ^ The LZSM problem is generally setup as 
follows. Define the eigenstates of a z to be |1) and \2) and 
set <p = so that j3 = b x and a — 2b z ; when \b z \ » l&^l, 
these states are approximate energy eigenstates. Nonzero 
b x produces an anti-crossing with an energy gap of 2b x 
(see Fig. [5]) which may be time-dependent. Now suppose 
that we drive b z through the anti-crossing, starting from 
some large negative value at t = up to a large positive 
value at t = T. Assuming that the system is initially 
prepared in state |1) at time t = 0, the probability P2(T) 
for the system to be in state |2) at time t = T is 

P2(T) = | M21 (T)| 2 = sin 2 [ X (T)]. (12) 

The fact that this depends only on x(T) demonstrates the 
suitability of the present formalism to the LZSM prob- 
lem. If we choose x such that x(0) = and x(T) = 0, 
then we achieve a perfect LZSM transition: the system is 
driven through the anti-crossing and returns to state |1) 
with probability 1. On the other hand, we may choose 
x{T) = 7r/2, in which case the system undergoes NARP 
and ends up in state |2) after being driven through the 
anti-crossing. Other values of x(T) lead to superpositions 
of |1) and 1 2). We may also consider LZSM interferome- 
try, where the system is driven through the anti-crossing 
periodically, and the resulting time-averaged probabili- 
ties to be in state |1) or |2) after many periods is again 
largely determined by x(T), as we will see. In choosing a 
x(t) for the LZSM problem, we must impose appropriate 
initial conditions. For simplicity, we focus on the case 
b z (0) — — oo, b z (T) — oo, for which we need x(0) < 0, 
x(T) > 0; the analysis can be extended straightforwardly 
to the case where b z is finite at t — 0, T. 

For constant b x , a simple example which satisfies these 
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FIG. 3: (Color online) Control field and NARP probability 
from Eqs. |13]) , (12| , © with (a),(b) %{T) = k/2 and b x T = 
tt/2, 1.6, 1.7, 1.8, 1.9, 2 and (c),(d) x(T) = tt/4 and b x T = 
tt/4, 0.8, 0.9, 1,1.1,1.2. 



boundary conditions and the QSL constraint is 



X 



b x t 



ib v T 



^t\ 



(13) 



where choosing < a < 16/(3T 2 ) ensures that the b z = 
a/2 from Eq. (J8]) remains finite in the interval t € (0, T), 
and r\ = — 1. This x saturates the QSL constraint when 
a = 0, implying that the corresponding b z will imple- 
ment near-QSL evolution for small a. To achieve a tar- 
get evolution characterized by x(^); sett — T in Eq. f|T3]) 
and solve for a: a(T) = 6[b x T - X (T)]/(b x T 3 ). Plugging 
this expression into Eqs. (TLH) and © yields a family of 
driving fields b z parametrized by T that will achieve the 
desired evolution for any value of x(T) e (0,7r/2]; some 
of these solutions along with the corresponding NARP 
probabilities are shown in Fig. |3l The restrictions on a 
become restrictions on T through the expression for a(T): 
x{T) < b x T < 9\(T); while the upper bound stems from 
the particular choice of Xi Eq. (Ti~3"l) . the lower bound is 
the familiar, universal QSL and gives rise to the step-like 
curves in Figs. [3^ and |3j:. These curves reveal that the 
desired LZSM transition is achieved as quickly as possi- 
ble by first driving b z to zero very rapidly, allowing the 
system to evolve under (3(t) for a time T < Tqsl, and 
finally driving b z quickly up to its final, large value. In 
addition to performing NARP, these near-QSL driving 
fields could be important for LZSM-based generation of 
entanglement in superconducting qubits^ where fideli- 
ties are often limited by short relaxation times. 

In the context of LZSM interferometry, the formalism 
of Eqs. (O,© yields an exact formula for LZSM inter- 
ference patterns. To show this, we begin by constructing 
a periodic driving field b z of period 2T in which the first 
half of one period is described by some choice of x(t) and 
the second half by X 2(t) — x(%T — t), corresponding to 
b z retracing its path. Using Eq. ©, it is simple to find 
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FIG. 4: (Color online) Time-averaged excited state proba- 
bility Vi from Eq. (|15[) with periodic driving field generated 
by Eq. (|13p with x(T) — n/2.1 as a function of the drive 
half-period T. 



the evolution after one full period: 

ttn(2T) = e 2 ^( T )cos(2 X (T)), 

u a i(2T) = -tsin(2x(T)), (14) 

where £0 = (£+ + £-)/2 and we have assumed that 
x(T) — —r]/3(T) for simplicity. From this expression, 
it is then straightforward to compute the time-averaged 
probability of being in state \2) (which is the excited 
state since each period begins and ends to the left of the 
anti-crossing) after many periods: 



P 2 = [2 + 2cot 2 (2 X (T))sin 2 (2£ (T))]- 



(15) 



Thus, we see that the present formalism readily produces 
a general, exact, analytic formula for V2 , whereas ana- 
lytic expressions for this important quantity typically re- 
quire several approximations. 5 It is clear that this func- 
tion takes values in the range [0, 1/2] and that V2 = 1/2 
for X (T) = tt/4, while V 2 = for X (T) = tt/2. This is 
to be expected since x(T) = tt/4 corresponds to a per- 
fect 50-50 "beam splitter", while xCO = 7r/2 ensures 
the system remains in the ground state after every sweep 
through the anti-crossing. In the context of charge qubits 
in which an electron inhabits a tunnel-coupled double 
quantum dot and states |1) and |2) correspond to the 
electron being in either the left dot or the right, the case 
Vi = can be interpreted as complete coherent destruc- 
tion of tunnelin g 45 i 46 since the electron is then fully lo- 
calized in one of the dots. For more generic values of 
x(^)> ^2 is modulated by the phase £o(T), which can 
be tuned independently of x(T) and which can produce 
interference fringes as control parameters are adjusted. 
This is illustrated in Fig.@]for the example from Eq. (fTS"]) 
with x(^) = tt/2.1. As the half-period drive time T is 
increased, an interference pattern emerges in which the 
peaks of the pattern sharpen and eventually disappear as 
x(T) approaches ir/2M. Interestingly, Fig. U reveals that 
there is a peak at the QSL time T = Tqsl — Tf/(2.1b x ); 
this is generally the case since at the QSL, |x| = so 
that £o(Tqsl) = 0. This leads to the surprising conclu- 
sion that if one wishes to trap the system in state |2) 
by choosing x(^) = 7r /2, then driving very close to the 
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QSL may not be ideal since small deviations away from 
x(T) = 7r/2 will produce the peak at T — Tqsl and, 
hence, large uncertainty in the state of the system. 

In conclusion, a general formalism for deriving exactly 
solvable time-dependent two-level quantum systems was 
presented. This formalism can vastly increase the num- 
ber of known exact solutions for physical Hamiltonians, 
as has been demonstrated with explicit, important ex- 



amples. These examples clearly show that this theory 
is a powerful tool in the design of control pulses both 
for quantum computation and for precise Landau-Zener 
interferometry near the quantum speed limit. 

I thank Michael Berry, Lev Bishop, Lukasz Cywihski, 
and Sophia Economou for helpful discussions. This work 
is supported by LPS-NSA and IARPA. 
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